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^ Abstract: The Deep Inelastic Scattering 2+1 jet cross-section is a useful 

Q\ observable for precision tests of QCD, e.g. measuring the strong coupling con- 

stant a s . A consistent analysis requires a good understanding of the theo- 
retical uncertainties and one of the fundamental ones in QCD is due to the 
renormalisation scheme and scale ambiguity. Different methods, which have 
been proposed to resolve the scale ambiguity, are applied to the 2+1 jet cross- 
^ section and the uncertainty is estimated. It is shown that the uncertainty can 

be made smaller by choosing the jet definition in a suitable way. 

in 1 Introduction 

o 



Oh 



Deep inelastic scattering (DIS) [1] has since long been a way of probing the nucleon struc- 
ture and studying the theory of strong interactions, Quantum Chromo Dynamics (QCD). 
The hadronic final state has been studied at fixed target lepton scattering experiments 



where effects such as high-p^ particle production [2] and the onset of jet production [3] 
have been observed, but the energy has not been high enough to have a substantial rate 
of clear multi-jet events. With the order of magnitude increased cms energy in the HERA 
electron-proton collider at DESY (Hamburg) the study of these phenomena enter a new 
era. 

In fact, multijet events have already been observed at HERA [4] and with increasing 
statistics they will be useful for precision test of QCD. In particular, the inclusive rate of 
2+1 jet events (where +1 denotes the jet emerging from the proton remnant) can be used 
to extract the gluon density xg(x) in the proton [5] or to measure the strong coupling 
constant a s and observe its 'running' [6], i.e. a s = a s (/i#) where [i 2 R is the renormalisation 
scale. An advantage at HERA is that this running can be observed in a single experiment 
to avoid systematic uncertainties that may enter when combining measurements from 
different experiments, as has been necessary so far (see e.g. [7]). 

In order to perform such analyses the theoretical cross-section has to be precisely de- 
fined and calculated with a proper understanding of its uncertainty. This requires the 
complete next-to-leading order (NLO) cross-section formalism, which also facilitates a 



1 



well-defined meaning of a s . As with all field theories, QCD has to be renormalised to get 
finite answers for the predictions of physical observables. To perform the renormalisation 
one has to define a renormalisation scheme which defines how much of the finite parts are 
to be kept at each order in the perturbative expansion of an observable and a renormalisa- 
tion scale which defines at which point or momentum scale the subtraction of the infinite 
parts should be made. From a physical point of view the result should be independent 
of the renormalisation scheme and scale, the so called renormalisation group invariance. 
This independence will, however, only hold for the complete perturbative series, whereas 
a truncation to any given order introduces a dependence on both the renormalisation 
scheme and scale (in next-to-leading order a change in the renormalisation scheme is 
equivalent to a change in the renormalisation scale). 

To solve this ambiguity at least three different methods for choosing the renormalisa- 
tion scale have been proposed, the BLM, FAC and PMS scales. 

BLM or Automatic Scale Fixing as the inventors Brodsky, Lepage, Mackenzie [8] call 
it. The renormalisation scale is chosen so that the second order terms with Nf 
dependence (iV/=number of active quark flavors) are absorbed into the running of 
a s . 

FAC Fastest Apparent Convergence as proposed by Grunberg [9] is defined by requiring 
all second order contributions to be absorbed into the running of a s so that the 
NLO result is equal to the LO result. 

PMS Principle of Minimum Sensitivity as proposed by Stevenson [10]. The renormalisa- 
tion scale is chosen so that the scale dependence is minimal, i.e. dR/dfi 2 ^ = where 
R is the physical observable under consideration. 

It is the purpose of this paper to apply these alternative scale choosing methods to DIS 
and, in particular, use them in order to asses the theoretical uncertainty represented by 
this free scale. This is necessary in order to estimate the theoretical error in, e.g., a 
measurement of a s or Kqcd- Previously, only simple estimates of the uncertainty has 
been made by varying the renormalisation scale with arbitrary factors around the simple 
choice of Q 2 which has no theoretical justification. Here, we perform a dedicated study 
where the theoretical uncertainty is more properly obtained by comparing the results 
obtained with these better motivated scale choices. 

In addition to the renormalisation ambiguities one has similar problems associated 
with the factorisation, which defines the way to absorb soft and collinear initial state 
singularities into the structure functions. This causes a corresponding scheme and scale 
(fip) dependence for the truncated perturbative series. The PMS method has also been 
applied to deal with these ambiguities, originally by Politzer [11] and later for example by 
Aurenche et al [12] who applied the PMS method to set the factorisation scale as well as 
the renormalisation scale for large p± processes involving real photons. In this dedicated 
study of the renormalisation scale no attempt will be made to optimise the factorisation 
scale, instead there will only be some comments on the applicability of the other methods 
to resolve that ambiguity. 

The paper is organised as follows. In section 2 we discuss the complete O(al) DIS 2+1 
jet cross-section and in section 3 we treat theoretical aspects of different renormalisation 
scale choices and comment on the factorisation scale dependence. The numerical results on 
scales and cross-sections are presented and discussed in section 4 and finally we summarise 
and conclude in section 5. 
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2 DIS 2+1 jet production 

2.1 Kinematic variables and jet definition 

The process of interest is (with four-momenta) 

electron(k) + proton(P) — > electron(k') + jetl(pi) + jet2(p 2 ) + remnant(p r ) 
as illustrated in Fig. 1. 
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Figure 1: Generic diagram for 2+1 jet production in DIS 
The overall ep cms energy squared is 

s = (P + k) 2 

and the inclusive deep inelastic scattering can be expressed in the variables 

Q 2 = -q 2 = —(k - k') 2 
W 2 = {P + qf 

giving the squared momentum transfer and invariant mass of the hadronic final state, or 
alternatively by the scaling variables 
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describe the relative energy sharing of the two jets from the hard scattering and 

is the incoming parton momentum fraction with respect to the proton (p = £P). In 
inclusive DIS, i.e. disregarding the hadronic final state, there are only two independent 
variables. These are usually taken as x,Q 2 or x,y and they are fully determined by the 
energy and angle of the scattered lepton. To describe the internal degrees of freedom in 
the 2+1 jet system, three additional variables are needed. These are usually taken as 
x p = Q 2 /2pq = z = z\ (giving z 2 = 1 — z) and the azimuthal angle <f> between the 
lepton scattering plane and the parton plane (p,p~i). 

Experimentally, jets can be defined according to the JADE- algorithm proceeding as 
follows. First the invariant mass of all pairs of particles is calculated. Then the pair 
with smallest invariant mass is identified and if its mass is below some cut-off mass the 
four momenta of these two particles are added to form a pseudoparticle (according to 
one of several possible schemes [7]). The process is repeated until all pairs of particles 
or pseudoparticles have invariant masses above the cut-off and the remaining particles or 
pseudoparticles are then identified as jets. The cut-off is given by 

{Pi +PjY = Sij > VcutW 2 (1) 

where y cut is a parameter between zero and unity that defines the jet resolution. Since the 
proton remnant (p r ) is counted as one of the particles in this procedure it is natural that 
the jet definition should scale with the total hadronic energy W 2 , but other choices are also 
possible. The proton remnant, which largely escapes in the beam pipe, can be represented 
by the momentum vector corresponding to the missing longitudinal momentum in the 
event . 

A direct comparison of experimental jets with hard parton emissions in matrix element 
calculations can be made if the same resolution criterion is applied in the calculation. 
Thus, the theoretical cross-section is formulated in terms of the invariant masses of pairs 
of partons. The target remnant, with momentum p r = (1 — £)P, is also included in 
order to automatically cut against collinear divergences in initial state emissions. The jet 
resolution y cut thereby defines the phase space for resolved parton emissions and sets the 
limits when integrating the matrix elements. 

2.2 Cross-section formalism 

Calculations of QCD matrix elements for jet production in DIS have been made since 
several years resulting in tree level diagrams up to order a 2 (see [13] and references 
therein). For our purposes the complete 0(a 2 s ) DIS 2+1 jet cross-section is needed, 
which have only recently been calculated (for one-photon exchange) by Brodkorb, Korner 
and Mirkes [13]. Prior to this result, the partial result obtained by Graudenz [14] by 
contracting the hadronic tensor with the metric tensor was available. The following results 
are mainly based on the complete calculation as implemented in the DISJET [15] program, 
but comparison is made with the earlier result as available in PRO JET [16]. Since both 
calculations are performed in the MS scheme [17] our results for the scales also apply in 
that scheme. 
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Assuming pure photon exchange and integrating over the azimuthal angle <f> between 
the parton plane and the lepton plane the cross-section can be written in the following 
general form, 

where (Ju+l (or F 2 ) represents the helicity cross-section for unpolarised photons and &l 
(or Fl) longitudinally polarised photons (the symbol U stands for unpolarised transverse). 
The longitudinal structure function Fl is a linear combination of the structure functions 
Fi and F 2 , Fl = F 2 — 2xFi which is zero in the Quark Parton Model (QPM) with spin 
1/2 pointlike constituents (Callan-Gross relation). In terms of the hadronic tensor 
we have 

f 1 Qx 2 p \ 

d(T U+L = I + -Q2~P»Pv J H ^ 

d*L = ^P.PvH^ 

where p is the momentum of the incoming parton, p = £P and x p = Q 2 /2pq = x/£. Other 
linear combinations are also possible like the one used by Graudenz [14], 



da 2ira 



era 



1 + (1 - y) 2 , 6(l-y) + y 2 
o a 9 + 7, a o 



(3) 



dxdy x 2 y 2 s 

where a g and a represents the 'metric' and 'non-metric' structure functions respectively, 

d(?g = -g^H^ 
da = ^P»P»H^ 

The two alternatives are related by a g = 2ctu + l — 3<7£ and a = &l- 



3 Renormalisation scale choices 



To make the dependence on the renormalisation scale fiR explicit, the leading order (LO) 
and next-to-leading order (NLO) 2+1 jet cross section can be written 



dx dQ 2 

d(?2+i 
dx dQ 2 



A 



A 



Qs(/4) 
2tt 

Qs(/4) 
2tt 



l+(BN f + C + b hi^ 2 R /Q 2 )) 



a 



2tt 



(4) 
(5) 



where A, B and C are functions of (x,Q 2 ,y cut , fip) after the dependences on the internal 
jet variables {x p ,z,<j)) have been integrated out. The scale dependence of the strong 
coupling constant is given by a renormalisation group equation, the so called /3-function, 
which to second order is given by 



d 



a 



d In [i 2 R 2-7T 
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(6) 
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where 



b 
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'-N r _ 



3 c 



-N f 
3 f 

lN c N f - C F N f 



2b 



and N c is the number of colours in QCD, Nf is the number of light fermions and Cf = 
(iV c 2 — 1)/2N C is the Casimir factor. The commonly used standard solution to Eq. (6) is 



a 



2tt & m(/4/A 2 ) 



^lnln^/A 2 ) 



Mn(/4/A 2 ) 



(7) 



where A is the QCD scale parameter which sets the normalisation of the strong coupling. 
It is convenient to use the dimensionless rescaled renormalisation scale 



/4 
Q 2 



(8) 



such that a direct comparison is made to the naive choice of Q 2 as the scale in a s . As an 
example the renormalisation scale dependence of the LO and the NLO expressions for the 
cross-section are shown in Fig. 2. Although the result displayed is for specific values of 
(x,Q 2 ,y cut ), it illustrates a general behaviour. With changing (x,Q 2 ,y cut ) the curves are 
displaced, but their functional dependence of the renormalisation scale is unchanged. As 
an illustration the three different scale choices (calculated as discussed below) according 
to BLM, FAC and PMS are also shown in Fig. 2. The 'geometric' interpretation of the 
FAC and PMS scales can easily be seen whereas for the BLM scale there is no such simple 
interpretation. 
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Figure 2: The renormalisation scale (p 2 = [i 2 R IQ 2 ) dependence of the 2+1 jet cross-section 
in leading order (dotted curve) and next-to-leading order (solid curve) at x, Q 2 and with 
y cut as indicated. The three different scale choices are also indicated. 
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3.1 The BLM scale 



The Brodsky-Lepage-Mackenzie (BLM) [8] method is inspired by QED where all fermion 
loops (vacuum polarisations) in the photon propagator are absorbed into the running 
of the coupling ct em . The same condition is applied to QCD after the renormalisation 
scheme has been chosen, i.e. all the quark and gluon vacuum polarisations in the gluon 
propagator are absorbed into the running of a s . In NLO it is sufficient to keep track of 
the Nf dependent terms since the vacuum polarisations have the form b K , so the iV c -part 
follows automatically with the iV/-part. This means that all Nf dependent terms showing 
up are a sign of bad bookkeeping and they should be absorbed into the running of a s . 
Although this principle is general it has to be applied with some caution in processes with 
gluon-gluon scattering in LO [8]. For the DIS 2+1 jet cross-section, there is however, no 
problem to find the renormalisation scale with the BLM method. The cross-section in 
Eq. (5) provides the condition 

BN f - l -NfH^ R IQ 2 ) = Q (9) 
from which the BLM scale is obtained as 

P 2 blm = exp(35) (10) 



Figure 3: The squared amplitude for the Nf dependent term in the NLO 2+1 jet cross- 
section. 

To O(al) the only process which gives an explicit iV/-dependence is illustrated in 
Fig. 3. The analytic form of this term can be obtained from the explicit cross-section 
formulas in [13, 14] and it is given by 

^ Jdx p dzd(/)A compton [| In (min{^-^(l - z),^y cut j^j - §] 

Jdxp dz d<j) (A com pt on -\- Af us i on ^ 
Where A compton is the Born cross-section for the QCD Compton process and 

A f us ion is the 

Born cross-section for the boson-gluon fusion process as seen in Fig. 4. (Comparing with 
Eq. (5) we have A = fdx p dzd(j>(A compton + A fusion ).) 



3.2 The FAC scale 

A common way of judging the reliability of a series expansion in physics is to look at the 
relative size of the last term calculated. If this term is small then the expansion is consid- 
ered as being trustworthy and if it is big then it is not. The method of Fastest Apparent 
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(a) 



(b) 





Figure 4: Illustration of the two Born processes for 2+1 jet production on parton level: 
(a) QCD Compton or gluon radiation, (b) Boson-gluon fusion. 



Convergence (FAC) [9] makes use of this criterion by choosing the renormalisation scale 
so that the last term in the perturbative expansion equals zero. In our case this means 
that the FAC scale is obtained from the requirement that the a 2 term in Eq. (5) is zero, 
i.e. 

BN f + C + b o \n( f i 2 R /Q 2 ) = (12) 

which gives 

2 _ / BN f ('■ 



Pfac ~ ex P 



bo 



(13) 



3.3 The PMS scale 

The Principle of Minimal Sensitivity (PMS) [10] defines the renormalisation scale by 
demanding that dR/ddfi 2 ^ = for any physical observable R. This choice is based on 
the knowledge that the complete perturbative series is renormalisation scale independent 
so one may try to impose the same condition on the truncated series. As showed by 
Stevenson in his original paper [10], demanding dR/dfi 2 ^ = will also make the optimised 
NLO result for R renormalisation scheme independent. In general the PMS principle 
should be applied both to the renormalisation scale and scheme. In that respect the PMS 
method differs from the BLM and FAC methods which only set the scale. 

A convenient way of parameterising different renormalisation schemes is through the 
coefficients (&;) in the /3-function, 

-f-^ = -b X 2 \l + WX + b 2 X 2 + ...] (14) 

where we have introduced the short hand notation 

A = (15) 

The first two coefficients, b and b\ are renormalisation scheme independent but the rest 
of them (bi,i > 1) are not. The PMS principle then amounts to requiring 

dR 

v = 

M = °" >1 
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The last condition is only applied for the b{ which are kept in the /3-function (the rest of 
the 6;'s are set to zero). In our case this means that only the first condition has to be 
applied to the cross-section. Taking the derivative of Eq. (5) and using Eq. (6) gives 



Pi„nlo 



d]n 2 -Ab \ 3 [W + 2(1 + M) (BN f + C + Mn(/4/Q 2 ))] (16) 



The PMS scale is then obtained from the requirement 



= ( 17 ) 

which gives 

W + 2(1 + WX) (BN f + C + b ln(/4/ Q 2 )) =0 (18) 

This is a transcendental equation so its solution can not be given as a closed expression 
but it can be solved numerically to arbitrary precision. For small A, the //^-dependence 
is mainly through the term b ln(/ifj/(5 2 ) in comparison to the b\X term, and one can 
therefore obtain the following iterative solution, 

A (n+1) = a '(Q 2ex p[- C ~ 260(1+6! A("))]) ^ 

with A(°) = 0. We have used the second iteration as the solution, i.e. 

2 / BN f + C W \ 

^ S = eXP ( bo 2fc(l + M«)] (20) 

which is precise enough for our purposes. Numerically we have (for Nf = 5) b = ^ 
and b\ = ||, so for the small A ~ 0.02 — 0.03 we are interested in the error is negligible. 
In passing we also note that comparing Eq. (13) and Eq. (20) we get the approximate 
relation p 2 PMS ~ 0.72p 2 FAC (for Nf = 5). 



3.4 Analytic expressions for the scales 

It would be instructive to have useful analytic expressions for the scales, but this is 
prohibited mainly by the long and cumbersome expression for C. We have therefore 
resorted to numerical calculation of the parameters A, B and C to get the different 
scales. 

There is also the question whether the scales should be calculated for the full 2+1 jet 
cross-section or the two Born processes separately. If one divides the cross-section into 
two parts, one for the qq final state and one for the qg, the following expressions for the 
BLM scales are obtained 

2 /jdx p dzd(/)A compton \n (min {^y cut , ^^(1 - z)}) $\ 

' comp ° n \ Jdx p dz d</> A compt 

on ' J I 

2 _ -1 

HBLM, fusion — L 

In the region of phase space where 1 ~ Xf (1 — z) does not contribute we get the interpre- 
tation that the renormalisation scale is simply given by the cut-off in the jet definition, 
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VcutW 2 , times a numerical factor. Theoretically one can have different scales for different 
processes. In fact, there is nothing that says that one should even have the same scale 
in the a s and a 2 contributions. In principle one could therefore also investigate the scale 
dependence of these two Born processes separately. We have chosen not to pursue this 
line mainly because experimentally one is first of all interested in the inclusive 2+1 jet 
cross-section, such that the extra uncertainties of quark and gluon jet identification does 
not enter. Furthermore, the use of different scales for the two Born processes will make it 
difficult to find a common y cut where the theoretical uncertainty is small. Having differ- 
ent y cut for the two Born processes would be very awkward experimentally, since it would 
require a separation of gluon jets from quark jets on an event-by-event basis before the 
jet reconstruction with resolution y cut has been made! 



3.5 The factorisation scale 



The dependence on the factorisation scale fip comes both from the parton densities and 
explicit a 2 terms of the form D\nfi F /Q 2 . Thus, the Oi^a 2 ) DIS 2+1 jet cross-section 
takes the general form 



dx dQ 2 
where 



1 + ( B{fi)N f + C(/4) + Mn(£|) + (D(fi.)bo + Etf F j) ln(^) ) A 



Q 



Jdx p dz d<f> A fusion (ii 2 F ) 



Q 2 



(22) 



(23) 



JdXp dz d(f> (A comp t on (flp^ + Af us i on (flp^ 

The factorisation scale dependence of the parton density functions is given by the GLAP 
[18] evolution equations 



d/ g (£,/4) 
d In fip 

din fip 



<x,(/40 f 1 

2"7T J(, Z 

<x,(/40 f 1 

2"7T J(, Z 



P qq (l) /,(*,a40 + p<g (A f G (z,fip) 



^ f q (z,fi 2 F ) + P GG f^)/ G (z,/4) 



where P qq etc are the splitting functions which are known to NLO [19]. With this in- 
formation it should be possible to minimise the factorisation scale dependence of the 
cross-section and therefore extract the PMS scale. However, for the FAC and BLM scales 
one needs some kind of analytic form for the structure functions so that the a 2 - and 
iV/-contributions can be identified. In this first step towards a complete understanding of 
the theoretical uncertainties in the DIS 2+1 jet cross-section we have limited ourselves to 
the renormalisation scale dependence. Only some preliminary numerical exercises will be 
made towards a study of the factorisation scheme and scale dependences that is needed 
to make the picture complete. 



4 Numerical results 

In order to calculate the renormalisation scales and the corresponding 2+1 jet cross- 
section the quantities A, B and C in Eq. (5) have been calculated numerically. This has 
been made for fixed (x,Q 2 ,y cut ) using the Monte Carlo program DISJET 1.0 [15], which 
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is based on the complete 0(ct 2 s ) calculation [13]. For comparison we have also used the 
Monte Carlo program PROJET 3.6 [16] which only includes the a 2 term proportional to 
cr g in Eq. (3). The calculations were made for the HERA design energy, 30 GeV electrons 
on 820 GeV protons giving a cms energy of y/s = 314 GeV. The dependence on the main 
DIS variables x,Q 2 and the jet resolution cut-off y cut will be of primary importance. We 
have therefore chosen (a:, (5 2 )-points to cover most of the kinematic range accessible at 
HERA as can be seen in Fig. 5. The ^-values considered are 0.0001, 0.0003, 0.001, 0.003, 
0.01, 0.03, 0.1, 0.3 and 0.9 and the (^-values are 9, 90, 900 and 9000 GeV 2 . 

10 5 

Q 2 

10 4 
10 3 
10 2 
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Figure 5: (x,Q 2 ) -points in the HERA kinematic range, with polar angle acceptance limits 
for the scattered electron and quark-jet (in QPM). 

The y cut values have been chosen to cover what can be considered as a possible region, 
i.e. y cut = 0.005,0.01,0.02,0.04,0.08. For y cut smaller than ~ 0.01 one expects the 
perturbative series to break down due to large logarithms of y cut causing a s In y cut <C 1 
to be invalid. This problem can in principle be solved by a resummation of such terms 
(similar to the leading InQ 2 resummation in the GLAP equations.) For y cut larger than 
~ 0.08 [15] the neglect of terms proportional to y cut in the cross-section calculations is no 
longer expected to be valid. 

For the required proton structure functions we have used parton density parameteri- 
sations from the package PAKPDF 2.1 [20], where we select the latest parameterisation 
by the CTEQ collaboration [21], set CTEQ2M which is a NLO fit including the first 
HERA data. The structure functions are given in the MS scheme and for the QCD scale 
parameter we have used their fitted value = 0.213 GeV . The choice of parton density 
functions should not play a big role for the resulting scales and cross-sections since the 
jet definition introduces a cut-off against small momentum fractions £, 

t > x + y cut (l - x) (24) 
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and all modern parameterisations more or less agree for £ > 0.01. The factorisation scale 
has been set to Q 2 which is arbitrary, but convenient for our purposes. Some estimates 
with other choices will be given below. 

4.1 Resulting renormalisation scales 

A selection of the scales obtained using the complete second order 2+1 jet cross-section 
(DISJET) is shown in Fig. 6. For a complete collection of all the scales calculated see 
Table 3 in appendix A. 

Fig. 6 demonstrates that the scales obtained with the different methods may differ 
quite substantially. The FAC and PMS scales are quite close, as expected from their 
definition (and discussed at the end of section 3.3). The BLM scale has, however, a 
different behaviour and in some (x,Q 2 ,y cut ) points it differs by orders of magnitude. The 
following trends are clear from studying Fig. 6. 

• The BLM scale decreases with increasing a;, whereas the FAC and PMS scale do 
not have any simple x dependence. 

• All three methods gives [i 2 R approximately proportional to Q 2 . 

• The BLM scale increases with increasing y cu ti whereas the FAC and PMS scales 
increase with increasing y cut for small x but increase with decreasing y cut for large 
x. 

Comparing the scales one sees that they cross around y cut = 0.01 for small x and around 
Vcut = 0.04 for large x. The increase of the FAC and PMS scales observed for large 
x at small y cut is unphysical and indicates that the perturbative expansion is breaking 
down. On the other hand, the BLM scale behaves in a way which is what we expect 
physically, it increases when y cut is increased and decreases when x is increased. This is 
a natural consequence of the jet-definition which is a measure of the typical virtualities 
in the process. We are therefore led to believe more in the BLM scale than the FAC and 
PMS scales even though the perturbative expansion can be questioned when the PMS 
and FAC scales start to grow in this way. 

The increase of the FAC and PMS scales at small x and large y cut is different from what 
has been obtained by Kramer and Lampe [22] for 3-jet rates in e + e~ at LEP energies. 
Their result seems to be more like what is obtained for large x ~ | as seen in Fig. 7. At first 
one might expect that the results in DIS should be similar to e + e~ since the processes are 
related on parton level through crossing. There are however some important differences. 
The e + e~ calculation is obtained by contracting the hadronic tensor with the metric one. 
This means that the 'non-metric' part defined as in Eq. (3) is not present in the e + e~ 
calculation. Furthermore, at small x the boson-gluon fusion process starts to dominate 
in DIS due to the increase of the gluon density. It is presumably the large contribution 
of this process to the 'non-metric' longitudinal part which is the main difference between 
DIS and e + e~ in this context. 

The scales have also been calculated using PROJET [16] which only includes the 
contraction of the hadronic tensor with the metric tensor for the 0(ct 2 s ) cross-section. 
Therefore only the part proportional to a g in Eq. (3) can be used to determine the scales. 
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Figure 7: Comparison of the scales obtained from the BLM (solid), FAC (dashed) and 
PMS (dotted) methods for (a) 2+1 jet production in DIS at x = 0.3, Q 2 = 9000 GeV 2 
and (b) 3-jet production in e + e~ as given in ref. [22]. 
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Figure 8: Comparison of the scales obtained in x = 0.001, Q 2 = 90GeV 2 according to the 
BLM (solid), FAC (dashed) and PMS (dotted) methods with (a) the complete a 2 cross- 
section (DISJET) and (b) the partial result (PROJET). 



The values for the BLM scale are almost the same as those obtained with DISJET but the 
FAC and PMS scales are quite different. As an example the scales obtained for x = 0.001 
and Q 2 = QQGeV 2 are shown in Fig. 8. The differences observed can be understood 
from the fact that the a 2 corrections to the 'non-metric' part (<r ) of the cross-section 
are relatively larger than the a 2 corrections to a g . The differences at large x are much 
smaller than at small x which is also natural since the 'non-metric' longitudinal part of 
the cross-section is relatively larger at small x and therefore gives larger impact there. 
For the BLM scale, on the other hand, the iV/-dependent terms are the same for the 
unpolarised and longitudinally polarised contributions, so that the resulting scale is the 
same. Comparing the scales obtained with PROJET with the ones from e + e~ one sees 
in this case a closer agreement for all a;, as expected when the 'non-metric' part of the 
cross-section is neglected. 
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4.2 Cross-sections with different scale choices 



With the scales calculated it is easy to get the cross-sections from Eq. (5) for the respective 
scales. The results are shown in Fig. 9 for the same (a:, (5 2 )-points as for the scales in 
Fig. 6. The cross-section is seen to increase with increasing x and decreasing y cut as a 
consequence of an increased phase space for resolved jets as given by the jet resolution 
(sij > Vcut^rQ 2 )- The strong decrease of the cross-section with Q 2 is mainly the trivial 
propagator effect from the photon exchange. 

In general we observe smaller differences between the cross-sections than between the 
scales and the decrease is more than just due to the logarithmic scale dependence of 
the running coupling constant. This is also what we expect from Fig. 2 where the scale 
dependence of the LO and NLO cross-sections can be compared. However, for large x 
and small y cut the cross-section for the BLM scale diverges as the a 2 corrections become 
negative and larger in magnitude than the Born term. 



4.3 Uncertainty due to renormalisation scale 

It is important to estimate the error in the theoretical calculation of the cross-section in 
order to use it for precision measurements of a s etc. There is, however, no 'correct' way of 
doing this and one is therefore left with some arbitrariness or prejudice (unless the exact 
answer is known). In our case we only know the first two terms in a series that may even 
not converge properly (see however the fourth paper listed in [10] for a discussion on the 
possible convergence of 'optimised' approximations and the Borel summability of QCD 
and also [23] for a comparison of the NLO and NNLO optimised scales in R e + e -). 

As already stated, the FAC method is based on one of the more common methods for 
estimating errors, the apparent convergence of the series. If this criterion is used then the 
error would be zero for the FAC scale, f or the PMS scale and (C + ^B)\ BLM 

for the BLM scale. One could also be very conservative and compare the three LO and 
NLO cross-sections simultaneously and assign the global difference as an error. This is, 
however, unnecessarily pessimistic since it ignores some knowledge. 

Another possibility is to look at the logarithmic derivative of the cross-section with 
respect to the renormalisation scale, din a/dln [i 2 R . By definition the error would then be 
zero for the PMS scale, b\\ 2 FAC for the FAC scale and 



9 In, 



d In ft 2 R 



= *blm [61 + 2(1 + WXblm) {C + fB))] 
BLM [i + (c + f b) \ blm ] 



for the BLM scale. If the renormalisation scale dependence is small, then the situation is 
at least more favourable than if it is large. In the latter case it is clear that the theoretical 
uncertainty is high, whereas in the former one it can be small even though this cannot be 
said with any certainty. In other words, it is a necessary but not sufficient condition that 
the renormalisation scale dependence is small. 

It has been proposed by Brodsky and Lu [24] that din a/dln ^ 2 r\ blm {i.e. evaluated 
at [i 2 R = PblmQ 2 ) 1S a suitable measure of the theoretical uncertainty. The argumentation 
for this is as follows. First of all the BLM scale should be chosen since this is the only 
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x = 0.0001 x = 0.001 x = 0.01 x = 0.1 



scale which behaves in a physical way. Secondly din a/dln [i 2 R is a measure of the renor- 
malisation scale dependence which indicates the size of the uncalculated terms. Therefore 
this has to be small and the ideal case would thus be when the BLM and PMS scales 
coincide. 

In Fig. 9 this error estimate is plotted for the points referring to the BLM scale and 
it is also given in Table 1. Both the figure and the table shows that the theoretical 
uncertainty becomes larger when Q 2 is decreased and x and y cut are held fixed. This 
cannot be concluded from just comparing the scales in Fig. 6. The explanation is mainly 
that the logarithmic derivative (Eq. (25)) is proportional to a 2 so that the running of the 
coupling gives this effect. 



dlntr/dln fi 2 R \ BLM 



X 


Q 2 

[GeV 2 ] 


Vcut 


.005 


.01 


.02 


.04 


.08 


0.0001 


9 


0.050 


0.032 


0.027 


0.185 


0.917 


0.0003 


9 


0.027 


0.033 


0.004 


0.110 


0.440 


0.001 


9 


0.038 


0.018 


0.016 


0.046 


0.215 


0.001 


90 


0.024 


0.003 


0.001 


0.029 


0.103 


0.003 


9 


0.238 


0.029 


0.020 


0.006 


0.101 


0.003 


90 


0.085 


0.014 


0.009 


0.002 


0.048 


0.01 


9 


3.027 


0.256 


0.021 


0.013 


0.018 


0.01 


90 


0.334 


0.065 


0.001 


0.013 


0.006 


0.01 


900 


0.132 


0.032 


0.001 


0.004 


0.009 


0.03 


9 






0.227 


0.003 


0.028 


0.03 


90 


4.257 


0.282 


0.046 


0.013 


0.013 


0.03 


900 


0.605 


0.114 


0.020 


0.007 


0.009 


0.1 


90 






0.366 


0.032 


0.022 


0.1 


900 




0.808 


0.114 


0.014 


0.012 


0.1 


9000 


55.21 


0.271 


0.068 


0.014 


0.001 


0.3 


900 






0.621 


0.069 


0.019 


0.3 


9000 




1.693 


0.170 


0.032 


0.012 


0.9 


9000 




1.442 


0.134 


0.174 


0.164 



Table 1: Renormalisation scale uncertainty as defined by din a / din fi R \ BLM , with bold 
face to show the region with better than 5% theoretical precision. (Missing values indicate 
negative cross-section, cf. Fig. 9.) 

4.4 Preferred ?/ cur values 

With the theoretical uncertainty defined in the above described way one can exploit its 
dependence on the jet resolution (y cu t)- It is easy to see from Table 1 that by choosing 
y cut in a suitable way one can always make the theoretical uncertainty small. Letting, 
e.g., <9 In (t/<9 In /4 1 

blm ^ define a small theoretical error gives the acceptable y cu t 
values indicated with bold face in Table 1, i.e. for small x one should use y cut ~ 0.01 and 
for large x one should use y cut ~ 0.04. 

One can also see from Table 1 that there is an interval of useful y cut values, especially 
at small x there is a clear upper limit which is a new feature previously not encountered. 
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The lower limit can be understood from the need to resum large logarithms of y cut and the 
upper limit is probably due to the neglect of terms proportional to y cut in the calculations 
which makes the a 2 corrections too much negative. 

We note that this dependence on y cut is actually beneficial from an experimental point 
of view. First of all, a smaller y cut gives a larger rate and thereby statistics. Secondly, 
to be able to extract the gluon density at small £ [5] one wants to use a y cut as small as 
possible in accordance with Eq. (24). 



4.5 BLM scale: parameterisation and analytic expression 

Since the BLM scale has been shown to be almost strictly proportional to Q 2 it is natural 
to perform a fit in a; with all (5 2 -points taken into account. The fit has been made for 
Vcut = 0.02 since this seems to be the best compromise if one insists on using the same 
y cut in the whole kinematic range. The following function has been fitted 

p 2 = exp(A + B\nx + C\n 2 x + D\n 3 x) (26) 

with the result A = -9.25, B = -3.37, C = -0.425, D = -0.0191. The fit can bee seen 
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Figure 10: Fit (solid) of p 2 BLM using Eq. (26) for y cut = 0.02 using all x , Q 2 -points together 
with the analytic expression Eq. (27) (dotted). 

in Fig. 10 together with the analytic expression 

2 ( 5\ 1 — x 

p = exp I -- I Vcut— — (27) 

The latter is expected to hold when the Compton process dominates and the term 1 ~ Xf (1 — 
z) does not contribute in the integration in Eq. (21). The deviations from these conditions 
explain the behaviour of our numerically calculated BLM scale. At small x the fusion 
diagram dominates and therefore we expect the BLM scale to be close to Q 2 . At large 
a:, the expression 1 ~ Xf (1 — z) is smaller than -^ry C ut in parts of the phase space and 
consequently the BLM scale gets smaller due to minimum condition in Eq. (21). 
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4.6 Comparison of Q 2 and the BLM scale 



To see the potential effects of using the BLM scale instead of Q 2 as renormalisation scale, 
the fit to the BLM scale given above has been used to calculate the 2+1 jet cross-section 
in two (a;, (5 2 )-areas for y cut = 0.02. The calculation has been made with the DISJET 
program for a high-Q 2 and a \ow-Q 2 bin with the following limits, 

low-Q 2 : .0001 < x < .02 
.05 < y < .95 
10 < Q 2 < 100 



high-Q 2 : .001 < x < .1 
.05 < y < .95 
100 < Q 2 < 9400 

The resulting cross-sections are given in Table 2 and show very small differences (~ 1%) 
between the two scales. 

o- 2+ i [pb] 





low Q 2 


high Q 2 


/4 = Q 2 


1269 


472 


f'R = PbLmQ 


1257 


479 



Table 2: Cross-sections, in pb, for the low and high Q 2 bins. 



However, it is important to keep in mind that even if o-(p, R = Q 2 ) ~ o-(p, R = PblmQ 2 )i 
the Kqcd extracted from analysing the running of the coupling will be different in the 
two cases. As an illustration, consider a measured 2+1-jet cross-section in one (x,Q 2 )- 
point. With the quantities A, B and C calculated it is then possible to extract a s once 
the renormalisation scale is chosen. To a first approximation the extracted value of a s is 
independent of [i 2 R since the sum BNf + C + b ln(/i|j/ Q 2 ) is only weakly dependent on [i 2 R . 
With a s given, the value of Kqcd will therefore scale with p since a s ~ 1/ hi(fi R / Aq CD ). 
To be more quantitative and taking the /^-dependence into account, assume that 

j^^(x = 0.1,Q 2 = 900,y cut = 0M) = 1782 pb/GeV 2 

has been measured. Depending on whether Q 2 or PblmQ 2 1S use( i as scale one gets the 
following results 

fi 2 R = Q 2 : a s = 0.1307 A^=213MeV 
/4 = P 2 blmQ 2 ■ oc s = 0.1370 AgL = 7 AMeV 

where p\lm = 0.05432 has been used. One sees that taking the /^-dependence from 
b \n(p, R / Q 2 ) into account makes Kblm > PblmKq2 = 50 MeV. 
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4.7 Factorisation scale effects 



Above we have only considered the = Q 2 . A complete analysis of the scale 

ambiguities in the DIS 2+1 jet cross-section would also have to include the effects of the 
factorisation scale uncertainty. Referring to Eq. (22) we see that fip ^ Q 2 will introduce 
extra terms, both iV/-dependent and iVf-independent, in addition to the /independence in 
the parton densities. A change of the factorisation scale will influence the renormalisation 
scale to an extent that depends on the factorisation scheme through the absorption into 
parton densities. As a rough estimate of the influence of the factorisation scale on our 
renormalisation scale results we have simply redone the determination of the latter by 
using Eq. (22) with fi 2 F = 0.1 Q 2 and fi 2 F = 10Q 2 . 

In short the effects can be described as follows. First of all the general (x,Q 2 ,y cut )- 
dependencies are not affected, but the renormalisation scales are shifted when the factori- 
sation scale is changed. For large x the effects are small but for small x they are quite 
considerable (factor 4-5). This big effect is more of a worst case than a realistic estimate. 
Physically we expect the factorisation scale to be of the same order as the renormalisation 
scale and since /tfj ~ Q 2 we do not expect the factorisation scale to be very different. As 
indicated, part of the /^-dependence should presumably be absorbed against other terms 
from the structure functions and thus not contribute or affect the renormalisation scale. 
It is clear, however, that this deserves a separate study to obtain a complete estimate of 
the theoretical uncertainty. 

5 Summary and conclusions 

In this paper we have investigated the renormalisation scale dependence of the DIS 2+1 
jet cross-section and three different methods for how to fix the renormalisation scale, the 
BLM, FAC and PMS methods. The different scales are obtained in the following ways: 
the BLM scale by absorbing all iV/-dependent terms into the running of a s , the FAC 
scale by requiring that the LO and NLO cross-sections coincide and the PMS scale by 
minimising the renormalisation scale dependence of the cross-section. 

The scales were calculated using the complete 0(a 2 s ) cross-section with the jets defined 
according to the JADE algorithm, S{j > y C utW 2 , for representative x, Q 2 and y cut values in 
the HERA kinematic range. The scales are approximately proportional to Q 2 whereas the 
x and y cut dependencies are more complicated. The BLM scale increases with decreasing 
x and increasing y cut which is what we expect physically since the typical invariant mass 
of pairs of partons goes as y cu tW 2 = y C utQ 2 ^r from the jet-definition. The FAC and 
PMS scales also increase with decreasing x for large y cu ti but for small y cut the trend is 
the opposite. There is also an unphysical increase of the FAC and PMS scales for large 
x when y cut is decreased. This indicates that the BLM scale is a better choice from a 
physical point of view. 

The scales obtained using the complete Oi^a 2 ) cross-section are also compared with 
the ones obtained using the partial cross-section from the contraction with the metric 
tensor and earlier results on 3-jet production in e + e~ [22]. For large x the scales more 
or less agree and behave in a similar way as a function of y cu ti but for small x the FAC 
and PMS scales from the complete DIS 2+1 jet cross-section behaves in a different way. 
The difference is mainly due to the dominance of boson-gluon fusion and, in particular, 
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its contribution to the 'non-metric' longitudinal part of the cross-section which is only 
present in the complete Oi^a 2 ) DIS cross-section. 

To estimate the theoretical uncertainty due to the renormalisation scale dependence 
we use the logarithmic derivative of the cross-section with respect to the renormalisation 
scale evaluated for the BLM scale. By choosing an appropriate value of y cut it is possible 
to make the uncertainty defined in this way small. The range of useful y cut values is 
Vcut ~ 0.01 — 0.04 with the smaller value for small x and the larger one for large x. Thus, 
a reliable theoretical result sets not only a lower limit on y cu ti but also an upper one. 

Our results on useful y cut values can be translated into other jet-definitions which also 
use invariant mass cut-offs, e.g. S{j > y' cut M 2 . The suitable values of y' cut are then simply 
given by rescaling of our results, y' cut = j^y C ut- 

Choosing y cut values as indicated above the differences between the cross-sections for 
the different scales is very small, ~ 1%. This can be compared with the difference between 
the LO and NLO cross-sections using Q 2 as scale which is in the order of ~ 10%. It is 
also interesting to note that for these y CU £-values the preferred scales are in the order of 
/ifj ~ Q 2 for small x and [i 2 R ~ 0.1 Q 2 for large x. 

A fit of the BLM scale (p 2 BLM ) as a function of x has been made for y cut = 0.02 using 
all Q 2 points. This fit is then used to calculate the cross-section for a low and a high Q 2 
region in the (a;, (5 2 )-plane. The difference in the cross-section obtained using the BLM 
scale and the simple Q 2 as scale is small, ~ 1%. Nevertheless it is important to use the 
correct scale if the running of a s is to be analysed or a proper value of Kqcd extracted. 

A complete analysis of the scale ambiguities in DIS would also have to include the 
factorisation scale dependence. The effects of using fip ^ Q 2 have been estimated crudely 
by calculating the renormalisation scales for two different choices, fi F = 0.1 Q 2 and fi F = 
10 Q 2 . The conclusions on suitable y CU £-values and cross-section differences are the same, 
but the numerical values of the scales are changed. 

In conclusion, we have demonstrated that with an appropriate jet definition it is 
possible to control the renormalisation scale uncertainty of the DIS 2+1-jet cross-section 
such that precision tests of QCD are facilitated. 
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helpful discussions and to Yu. Dokshitzer and G. Kramer for comments on the manuscript. 
We would also like to thank E. Mirkes for helpful communications regarding the DISJET 
Monte Carlo program. 
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A Summary of renormalisation scales obtained 
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Table 3: Obtained renormalisation scales, p 2 = fi^/Q 2 . 
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